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An automated pipeline for asteroseismology based on the 
autocorrelation of stellar time series 
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ABSTRACT 

The autocorrelation of an asteroseismic time series has been identified as a powerful tool capable of providing measure- 
ments of the large frequency separations. The performance of this method has been assessed and quantified by Mosser 
& Appourchaux (2009). We propose now an automated pipeline based on it and describe its performance. 
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1. Introduction 

The mean large separation in a seismic spectrum is a key 
parameter in asteroseismology. Observationally, identifying 
the large separation proves the detection of solar-like os- 
cillation. It was the case with the first determination of 
the exact large separation of Procyon (IMosser et al. 1 998) 
even if the energy excess of the oscillation was not iden- 
tified, or of the CoRoT target HD 1 75726 for which n o 
mode identification has been possible (jMosser et al. 2 009). 
The measurement of the large separation opens then the 
determination of other parameters and the possibility of a 
large scientific output (stellar modeling, mass and radius 
measurement . . . ) . 

With the event of large set of data produced by CoRoT 
and Kepler, it appeared necessary to develop efficient tools 
able to determine the parameters that describe an os- 
cillation spectrum (jMathur et al. 20091 IHekker et al. 20091 
IHuber et al. 2009). All these tools search first in the power 
spectrum for an energy excess due to the solar-like oscilla- 
tions. Then, when this excess has been identified, the sig- 
nature of a solar-like oscillation spectrum is searched with 
different methods, all dealing with the detection of a regular 
pattern a la Tassoul (Tass ouT 1980|) . 

We have chosen a different way to perform the analysis 
of seismic time series. Instead of searching for the regular- 
ity of the oscillation pattern in the Fourier spectrum, we 
search for its physical cause. The comb structure of the 
solar-like oscillation pattern being due to the fact that all 
pressure modes propagate throughout the stellar diameter 
at the same phase velocity, we search for this signature 
in the autocorrelation of the time series. The autocorrela- 
tion correlates in fact any wave with itself after a double 
travel across the stellar diameter. Instead of computing the 
autocorrelation, we deal with the calculation of the Fourier 
spectrum of the Fourier spectrum. In order to perform an ef- 
ficient analysis, we comput e narrow window filtered Fourier 
spectra, as proposed by |Roxburgh fc Vorontsov (2006) 



Mosser & Appourchaux (2009) hereafter MA09, have ap 



plied and quantified the method, and presented its imme- 
diate avantage: the signature of the large separation occurs 
at a time delay much smaller than the mode lifetime, so 
that the autocorrelation of the time series is efficient since 
it stacks cophased signals. 

In this work, we address specifically the description of 
the pipeline based on the autocorrelation of the time series. 
We recall briefly in Section [5] how the envelope autocorre- 
lation function is calculated (the full description is given 
in MA09). Section [3] presents the main steps of the au- 
tomated pipeline for the determination of the mean large 
separation, then the determination of the mode envelope. It 
also shows how to measure the variation of the large sepa- 
ration with frequency. A few characteristics and properties 
are discussed in Sect. [4] 



2. The envelope autocorrelation function 

The envelope autocorrelation function (EACF) presented 
in MA09 is first normalized at 1 at null delay, then rescaled 
to the noise contribution. This rescaling is a crucial step 
for the automated analysis, since it gives the possibility to 
test the performance of the method independently of any 
modeling. A statistical test as the null hypothesis allows us 
to define the reliability of the detection. This rescaling also 
shows that a tiny peak in the EACF may be in fact fully 
reliable when compared to the tiniest signature of a white 
noise. 

With this scaling, and taking into account properties 
of solar-like oscillation spectra, MA09 have shown that an 
EACF greater than 8 corresponds to the rejection of the 
null hypothesis at the level 1 %. White noise being then 
unable to explain the autocorrelation signal, one can derive 
the positive detection of a signal. After correction of pos- 
sible artefacts, as the signature of the orbital frequency of 
CoRoT ( Auvcrgn e et al. 2009[ ) , this signal can be identified 
to solar-like oscillations. 
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3. Automated pipeline 

3.1. Determination of the mean large separation 

The mean large separation (Av) is denned as the mean 
value of Av(v), in a frequency range centered on v max as 
large as the mode envelope. This asteroseismic parameter 
is derived first. This means that it is not necessary to first 
estimate the background due to stellar activity, what makes 
the asteroseismic analysis direct and simple. 

In order to perform efficiently the search, scaling re- 
lations are presupposed between ^ max , (Av) and the full- 
width at half-maximum Sv cnv of the mode envelope. The 
relation between i/ max and (Av) follows the monomial law 
reported in CoRoT data (jHekker et al. 20091 MA09). In or- 
der to examine possible outliers, ±30 % variations around 
this law are permitted. We have checked that the relation 
between z/ max and (Av) is satisfied for all targets with an 
agreement better than 20%. The relation between Sv onv 
and (Av) is also derived from CoRoT data (MA09): a ratio 
Sv env /(Av) of 4 is first supposed for red giants, and 10 for 
solar-like stars. 

The systematic blind analysis consists in searching (Av) 
in a range around an initial guess value Av g . One computes 
the Fourier spectrum of the windowed Fourier spectrum 
with a Hanning filter centered on f max (A^ g ) and a width 
Sv env (Av g ). This gives the (possible) signature Av close to 
Av g , if validated by an EACF greater than 8 as presented in 
Sect. [U Large separations in the range from 0.5 to 250 /iHz 
can be exhaustively tested in 18 steps, with the guess values 
Av g in progression with a geometric ratio \/2. 

The mean value (A^) being firmly established, a more 
detailed analysis around the interesting frequency range 
gives then ^ max and Sv cnv . Nevertheless, at this stage, these 
parameters v mllx and Sv cnv are less precisely estimated than 
(Av). In fact, the value of i/ max does not strictly correspond 
to the maximum oscillation signal, but to the maximum 
EACF signal. The difference derives from the fact that the 
modes lifetime varies with frequency: a mode with low am- 
plitude but large lifetime may have a larger autocorrela- 
tion signature than a mode with a larger amplitude and a 
shorter lifetime. 



3.2. Mode envelope 

The precise determination of v max , Sv cnv and of the bolo- 
mctric amplitude of radial modes can now be performed in 
a spectrum smoothed with a cosine window of width (Av) . 
Then, the background power due to the stellar activity is es- 
timated, in a similar way as performed by other pipelines, in 
order to estimate the power excess. We note that the prior 
determinations of ^ max and 5v env make their refined deter- 
mination easy. The bolometric of ra dial modes is performe d 
according to the recipe described in Michel et al. (2009) 



3.3. Variation of the large separation 

If the EACF is high enough, the determination of the vari- 
ations of the large separation with frequency is possible 
without any mode fitting, with narrow filters. The exam- 
ination of the derivative dAv/dn of the large separation 
with respect to the radial order n is best performed with 
a filter width approximately equal to Sv cnv /2. In order 
to fully investigate the variation of Av(v), a filter width 



equal to 2 (Av) is well suited. MA09 have shown that a 
filter width equal to 0.75 (Av) allows the identification of 
the degree of a ridge, what can be useful for F targets 
dAppourchaux et~a l. 2008). The threshold level for analyz- 
ing Av(v) is lower than for the blind analysis, since it ben- 
efits from the prior measurement of (Av) and from the fact 
that Av(v) is now searched in a narrow range: the 1 % re- 
jection level is at 4.6 when searching for variation around 
±20% of (Av) with a 2-(Av) wide filter. 



4. Discussion 

Compared to other pipelines, the pipeline based on the en- 
velope autocorrelation function shows many qualities: 

- It is rapid; the complete operation takes less than 10 s 
for the blind analysis of a CoRoT time series with 4 10 5 
points, corresponding to a 150-day run at a 32-s sampling 
time. For the first Kepler 44-day long time series with a 
sampling about 30min, the blind analysis lasts 0.3 s per 
target. Most of these operations are direct calculations or 
fast Fourier transforms (fft). The use of fft simply requires 
the scaling presented in Section [5] to be performed at the 
exact frequency resolution. 

- It includes the determination of a reliable threshold 
level. Contrary to other pipelines, the performance and the 
determination of error bars do not depend on any scal- 
ing based on the analysis of theoretical models. Similarly 
to other pipelines, the EACF just presupposes a solar-like 
regular pattern. 

- It measures the mean large separation with a reliable 
and reduced error bar. For a solar-like star (resp. a red 
giant), the determination of (Av) at the detection limit is 
achieved with a relative precision of 0.7% (resp. 1.9%). 
This precision scales with the inverse of the EACF (values 
up to 600 have been measured with CoRoT). Since the large 
separation may significantly vary along the spectrum, the 
definition and the significance of the mean large may be 
not accurate enough to fully describe the rich output of 
this asteroseismic observable (see next point). In the same 
way, the definition of i/ max presupposes a Gaussian mode 
envelope, what appears to be not valid in all cases, due to 
the stochastic nature of the mode excitation. 

- It gives access to the variation of the large separation 
without any mode fitting. 
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